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i. 


The  Technical  Reports,  Numbers  3,  4  and  5  constitute  the  major  part 
of  the  final  report  for  project  N00014-67-A-0385-0004,  ONR  032-509,  since 
they  contain  the  findings  that  have  been  achieved  with  the  finite  element 
method.  In  summary  these  are: 

STRESS  RAISER 

The  anisotropic  elastic  properties  of  a  material  can  act  as  a 
stress  raiser.  The  elastic  constants  used  for  study  were  derived  from 
the  measured  elastic  compliances  for  alumina.  In  other  materials  the 
degree  of  anisotropy  may  be  larger  and  hence  more  important  si  tee  alumina 
is  not  the  most  anisotropic  material.  The  direction  of  maximum  aniso¬ 
tropy  for  alumina  is,  in  terms  of  angles  defined  in  Reports  3  and  4, 
a  =  0,  $  =  55°. 

NOTCH  FACTOR 

The  notch  factor  as  calculated  by  Neuber  for  hyperbolic  cross 
section  specimens  does  not  apply  when  the  material  is  anisotropic.  In¬ 
stead  of  1.45,  the  FEM  value  is  1.82.  This  latter  value  may  not  be 
correct  since  there  are  some  obvious  deficiencies  in  the  model  used.  It 
may  be  even  highe.-  since  the  stress  gradients  were  greater  than  the  model 
was  designed  to  handle. 

NON-UNIFORM  STRESS  STATE 

The  stress  increase  for  a  bicrystal  in  tension  was  found  to  yield 
stress  concentrations  at  the  boundary  as  high  as  1.5.  It  is  thus  seen 
that  the  assumption  of  uniform  stress  states  in  bicrystal  mechanical 
response  is  not  generally  accurate. 

MODEL  FOR  CERAMIC  BODIES 

A  model  for  the  study  of  the  mechanical  behavior  of  a  polycrystalline 
ceramic  body  has  been  proposed.  For  certain  selected  orientations  the 
model  generates  results  that  are  consonate  with  reality  and  shows  that 
the  position  of  maximum  shear  can  be  shifted  by  a  change  in  anisotropy. 
There  seems  to  be  a  limit  to  the  effect  anisotropy  can  make  on  a  system 
after  the  degree  of  anisotropy  has  passed  a  certain  value. 


When  the  nodel  was  subjected  to  the  case  of  different  orientation 
for  each  crystal  the  stress  gradients  which  resulted  were  higher  than 
could  be  adequately  studied  by  the  model.  A  further  modification  of 
the  model  would  be  necessary  for  the  proper  characterization  of  a 
ceramic  body  by  this  method . 

A  further  technical  report  will  be  submitted  when  the  measurements  of 
the  strain  distribution  in  an  alumina  bicrystal  have  been  completed.  Thus 
far,  the  "d"  spacing  for  the  (054)  diffraction  peak  has  been  measured  at 
systematically  selected  points  over  the  surface  of  the  bicrystal.  The 
residual  strain  has  in  this  way  been  partially  measured.  Similar  data  are 
being  obtained  for  the  (330)  diffraction  peak.  Preliminary  measurements 
with  an  applied  tensile  stress  have  shown  a  linear  relationship  between 
applied  stress  and  the  inter planar  spacing.  Thus,  the  apparatus  is  seen 
to  be  in  good  operating  condition.  With  additional  time,  the  measured  stra 
in  a  bicrystai  of  alumina  under  tensile  loading  should  be  completed.  A 
request  for  a  no-cost  extension  has  been  made. 
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A.  Introduction 

Ceramics  are  rapidly  increasing  their  importance  as  technological 
materials  in  a  number  of  important  applications.  Due  to  this  increased 
interest  and  also  because  of  the  commercial  availability  of  pure  metal-oxide 
powders  and  single  crystals,  a  great  deal  of  new  knowledge  of  their 
mechanical  behavior  is  being  collected.  However,  the  progress  in  under¬ 
standing  some  characteristics  of  the  mechanical  response  of  polycrystalline 
aggregates  of  metal  oxides  is  seriously  hampered  by  the  inherent  analytical 
difficulty  encountered  in  modeling  such  systems. 

In  recent  years,  the  problem  was  approached  from  the  bi-crystal  point 
(1  2  3) 

of  view  *  ’  in  order  to  gain  insight  into  the  process  of  brittle  fracture 
at  the  grain  boundaries.  This  form  of  failure  is  sometimes  observed  in 
polycrystalline  metal  oxides. 

Brittle  fracture  may  be  induced  in  suitably  shaped  tensile  specimens 
provided  that  a  sufficiently  high  stress  level  is  maintained  in  the  gauge 
section  of  the  specimen. 

A  tensile  specimen  generated  either  by  a  hyperboloid  of  revolution  or 
by  shaping  the  two  opposing  sides  by  hyperbolas  as  shown  in  Fig.  1,  will 
display  a  substantial  tensile  stress  at  the  focal  plane  while  maintaining 
sufficiently  low  stresses  at  the  gripping  ends. 

The  hyperbolic  tensile  specimen  shape  was  used  in  measuring  the 
intrinsic  brittle  strength  of  single  and  bi-crystal  MgO  ^ . 


2. 


a)  Axisymmetric  b)  Planar 

Fig.  1.  Hyperbolic  Tensile  Specimens. 

The  stress  distribution  in  hyperbolic  tensile  specimens  was  studied 
by  Neuber  ^  and  by  Leven  ^ .  A  stress  concentration  factor  or  notch 
factor  for  homogeneous  and  isotropic  specimens  was  calculated  from  an  exact 

theory  and  was  experimentally  confirmed  to  be  accurate  within  a  few  percent. 

(6) 

Some  recent  investigations v  however,  considered  the  influence  of 
elastic  mismatch  across  a  planar  boundary  in  rectangular  bi-crystal  tensile 
specimens.  It  ./as  shown  that  the  anisotropy  of  elastic  properties  as  well 
as  the  elastic  mismatch  produce  similar  effects  to  that  of  notches  in  the 
tensile  specimens. 

Various  attempts  have  been  made  ^  to  find  an  exact  and  general  solution 
for  the  problem  of  stress  distribution  in  anisotropic  bodies  bounded  by 
cylindrical  surfaces.  The  theoretical  approach  was  successful  only  for 
infinitely  large  bodies  bounded  by  planes  or  parabolic  surfaces.  The  intro¬ 
duction  of  the  necessary  finite  boundary  conditions  to  delimit  a  finite 
body  would  prevent  solution. 

An  attempt  is  made  here  to  formulate  the  problem  in  such  terms  that  an 
approximate  method  of  analysis  could  be  applied,  which  as  a  limit,  would 
converge  to  the  exact  solution. 


3. 


B,  Method  of  Analysis 


Since  a  closed  fora  exact  solution  is  not  available  for  anisotropic 
finite  bodies  at  the  present  time,  a  solution  based  on  the  calculated  stiff¬ 
ness  of  the  anisotropic  finite  continuum  was  sought.  This  method  of  analysis 
has  been  used  in  the  past  in  various  forms  such  as  "Moment  Distribution"  or 
"Deflection  Analysis"  of  structural  elements.  The  Finite  Element  Method 
(F.E.M.)  of  analysis  and  its  terminology  was  first  formulated  by  M.  J. 

Turner  et  al.  and  it  was  first  applied  in  its  present  form  to  an  elastic 
continuum  by  R.  W.  Clough  It  has  been  shown  repeatedly  that, 

while  the  method  of  analysis  is  approximate,  it  converges  to  the  exact 
solution  as  the  limiting  value. 

The  first  step  of  applying  the  F.E.M.  analysis  is  to  divide  the  finite 
elastic  continuum  into  polygonal  elements.  The  elements  are  usually  tri¬ 
angular  or  quadrilateral  in  shape.  The  elements  are  considered  to  be 
jointed  at  their  corners  (the  nodal  points)  and  the  stiffnesses  of  the 
individual  elements  are  assumed  to  interact  through  these  mutual  nod?I  points. 

The  choice  of  element  size  requires  some  experience  of  the  user  so  that 

the  possible  error  in  the  calculations  is  minimized  for  a  given  number  of 

elements.  The  model  used  is  shown  in  Fig.  2, 

(12) 

It  was  shown  by  Wilson  that  the  error  E  for  the  displacements  u^ 

of  a  nodal  point  j  is: 
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~2 

0  3  u 

2x.y. 

3  3  3-3 


x  y 


(1) 


where  xvy  are  the  coordinates  of  the  nodal  point  j. 
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The  error  term  demonstrates  the  influence  of  both  the  element  size  and 
that  of  the  rate  of  change  of  strain.  Since  the  rate  of  change  of  strain  is 
determined  by  the  geometry  of  the  problem,  the  proper  choice  of  element 
size  will  improve  accuracy. 

While  it  is  generally  true  that  convergence  of  finite  element  solution 
increases  with  decreasing  element  size,  the  gain  is  sometimes  offset  Dy  the 
increase  in  the  "round-off"  error  resulting  from  the  limited  word  length  of 
electronic  digital  computers. 

In  the  stiffness  formulation  of  F.E.M.  analysis  the  unknowns  of  the 
continuum  are  the  nodal  displacements.  The  basic  equation  which  relates 
nodal  displacements  to  the  nodal  forces  F  will  always  take  the  form 
in  matrix  notation: 

{F}  =  [K]  {u}  (2) 

where  [K]  is  the  stiffness  matrix  of  the  continuum.  The  global  stiffness 
matrices  are  defined  as: 

[K]  =  [B]^[D] [B]  x  volume  (3) 

In  Equation  (3),  [B]  and  [D]  represent  the  displacement-strain  and 
stress-strain  matrices  respectively. 

The  global  stiffness  matrix  is  symmetric  and  only  the  boundary  con¬ 
ditions  render  it  non-singular.  Since  the  stiffness  of  a  particular  nodal 
point  is  defined  by  the  combined  effects  of  only  those  elements  which 
contain  that  particular  nodal  point,  the  stiffness  matrix  is  large  and  is 
very  sparsely  populated.  An  efficient  computer  storage  scheme  is  a 
necessity  to  handle  large  problems  successfully. 
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Once  the  global  stiffness  matrix  is  developed,  the  unknown  nodal  dis¬ 
placements  are  calculated  from  equation  (2). 

The  element  strains  and  stresses  are  computed  from  the  relations: 


[ej  =  [B] [u] 

(4) 

[o]  =  [D][e] 

(5) 

The  continuity  between  the  elements  is  assumed  by  a  suitable  linear 
polynomial  expression  for  [B]  to  relate  displacements  within  the  tri¬ 
angular  elements  to  that  of  stress. 

The  stress-strain  relationship  in  terms  of  the  principal  technical  con¬ 
stants  E,  G  and  U,  (Young’s  modulus,  shear  modulus  and  Poisson's  ratio 
respectively)  is  given  in  general  form  (6). 
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(6) 


Y  =  -±-  a 

xy  G12  xy 

Isotropy  is  re-established  by  the  proper  choice  of  the  principal  technical 
constants,  i.e.,  =  E;  V12=V21  and  G12  = 

An  excellent  description  of  this  technique  is  given  by  Zienkiewicz  and 
Chung  in  their  book  .  The  actual  matrix  manipulations  for  any  sizable 
finite  body  can  only  be  accomplished  b>  digital  computers  of  storage  capacity 
in  excess  of  100  Kbytes. 
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The  computer  program  used  in  this  investigation  was  written  in  its 

(13) 

original  form  by  H.  D.  Dahl  and  was  extensively  modified  for  the 

present  purpose.  The  program  is  written  in  FORTRAN  IV  language  and  con¬ 
sists  of  a  main  routine  which  in  turn  calls  for  subroutines  during  the 
processing  of  input  data.  The  method  of  solution  of  the  large  number  of 
simultaneous  linear  equations  of  expression  (2)  is  an  iterative  technique 
based  on  the  Gauss-Seidel  over-relaxation  method.  The  computer  output 
consists  of  the  re-print  of  input  data,  calculated  nodal  point  displacements 
and  of  the  stresses  at  the  centroid  of  each  element. 


C.  Stress  Distribution  in  Isotropic  Specimens 

The  nature  of  stress  concentration  in  hyperbolic  specimens  is  now  well 

understood  due  to  the  pioneering  work  of  Neuber  in  the  early  1930’s.  It  was 
(4) 

shown  that  in  homogeneous  and  isotropic  flat  specimens  under  pure 
tension  the  principal  stresses  O u  and  a ^  may  be  calculated  readily  from 
the  equations: 


a 


a 


A  2  2 

=  cosh  u  cosv  (2  +  cos  v„-  cos  v) 

U  h  h2 

=  — 7*  cosh  u  cosv  (cos  V  -  COS^  ) 
v  ,  4  o 

h 


where  A  =  P 


sinvD 


v  +  sinv  cosv 
o  oo 


2  2  2 
and  h  =  sinh  u  +  cos  v 


(7) 


u  and  v  are  elliptic  coordinates,  related  to  the  cartesian  coordinates  by 
x  =sinhu  •  cosv 
y  =coshu  •  sinv 


(8) 


The  curves  with  u  =  const,  are  ellipses,  those  with  v  =  const,  are  hyperbolas. 
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The  surface  of  the  hyperbolic  specimen  is  given  by  v  =  1  v  .  The  waist  of 
the  specimen  is  defined  by 

x  =  u  =  o 


yv=v0  =  sinv0  =  a 

Equation  (7)  gives  the  principal  stresses  in  terms  of  p,  the  average  axial 
stress  in  the  focal  plane  of  the  specimen.  It  is  obvious  that  both  stresses 
have  their  maximum  in  the  focal  plane;  the  values  follow  from  Equation  (7)  by 
setting  u  =  0.  For  axisymmetrical  specimens  generated  by  a  hyperboloid  of 
revolution  the  equations  are  slightly  more  complex  and  for  the  focal  plane 
the  stress  equations  can  be  written  in  the  form: 


(0  )  =  — x  {Bcosv  +  C  [l-(a  +  2)]  cosv} 

u  u=o  ,2  1 

h 

+  cosv  +  Ccog2v  _  Aj 

4 

h 


(CO  =  ^  {Ceos  (a  -  1)  -  A  -v°?— --  } 
v  u-o  ,  z  1  +  cosv 

n 

+  {A  -  B  -  C  cos2v} 

h 

(Va=o '  72 lA  iTs;;  - B  cosv  +  c  cosv  <a  - 2! 

h 


2  2  2 
h  =  sinh  u  +  cos  v 


a  =  2  Cl  -V) 


A  =  C(a-l)  (1  +  cosv  ) 

0 


B  -  A  -  C  cos  v 


c  =  4(- 


1  ,  V 

1  +  cos  o 


V  =  Poisson’s  ratio 


1  +  (2-a)  cosv  +  cosv^ 
o  o 


The  elliptic  coordinates  used  here  are  defined  by: 
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x  =  sinlm  cosv 
y  =  coshu  sinv  cosw 
z  =  coshu  sinv  sinv 

The  surfaces  with  u  =  constant  are  ellipsoids, 

those  with  v  =  constant  are  hyperboloids 

and  those  with  w  =  constant  are  planes  through  the  x-axis. 

The  specimen  is  bounded  by  the  hyperboloid  v  =  vq. 

The  photoelastic  measurements  of  Leven  (5)  showed  that  the  numerical 
values  calculated  from  the  above  equations  are  within  a  few  percent  of  those 
measured  experimentally.  In  Fig.  3a  and  3b,  the  axial  and  transverse  stress 
distribution  are  plotted  in  terms  of  p,  the  maximum  nominal  stress.  The 

parameters  of  the  generating  hyperbola  for  both  flat  and  axisymmetrical 

?  2 

specimens  were  chosen  as  0.00098  for  a~  and  0.00124  for  b  respectively, 
where  a~  and  are  the  constants  in  the  equation  defining  the  surface  of  the 
flat  specimen  (hyperbola  symmetrical  to  the  y-axis) 


(ID 


Using  Neuber's  parameter,  vQ,  defining  the  same  surface,  one  can  write 
this  equation 

fair-)2  -fc^>2  -1  (12) 


sanv„  cosv 
o  o 


O  O 

Since  sin  vQ  +  cos  vfl  =  1  (13) 

one  has  .o  multiply  a^  and  b^  with  a  normalization  factor,  f^,  so  that  the 

2  2 
last  equation  is  satisfied,  (f  =  450,  f  =  21.5)  It  follows  that  sin  v 

o 

2 

=  0.44  and  cos  vQ  =  0.56.  Inserting  vQ  in  Equation  (7)  gives  the  curve  plotted 
in  Fig,  3a  as  Nueber's  exact  solution. 

The  procedure  for  the  axisymmetrical  case,  plotted  in  Fig.  3b,  is  analog. 
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The  above  analyses  were  repeated  using  the  "Finite  Element"  method.  For  the 
purpose  of  analysis  a  model  was  constructed  consisting  of  288  elements  and 
175  nodal  points  to  represent  the  transverse  half  of  the  tensile  specimen  as 
shown  in  Fig.  2.  The  numerical  results  of  the  F.E.M.  analysis  are  also 
shown  in  Figs.  3a  and  3b.  It  can  be  seen  that  the  "Finite  Element"  analysis 
produces  a  very  close  set  of  results,  within  about  5  percent  of  the  "exact" 
solution  of  Neuber.  Some  part  of  this  error  nay  be  due  to  the  fact  that  the 
triangular  elements  are  finite  in  size  and  they  are  located  slightly  on  one 
side  of  the  focal  plane.  The  other  possible  source  of  error  is  that  the 
strain  is  considered  constant  in  any  one  element.  Since  tne  rate  of  change 
of  strain  is  large  in  the  immediate  vicinity  of  the  focal  section,  even  a 
small  finite  element  size  introduces  some  error  in  the  actual  numerical  value. 
The  accuracy  may  be  increased,  especially  in  the  area  near  the  edge  of  the 
specimen,  by  refining  the  triangular  mesh  used. 

In  this  instance,  however,  the  stress  distribution  follonrs  a  continuous 
and  smooth  curve,  thus  a  relatively  accurate  extrapolation  of  F.E.M.  data 
points  to  the  hyperbolic  boundary  is  possible. 

The  "exact"  and  F.E.M.  analyses  of  isotropic  tensile  specimens  served 
usefully  for  checking  the  relative  accuracy  of  the  crystal  model  and  that  of  the 
computing  process. 

D.  Stress  Distribution  in  Anisotropic  Single  Crystal  Tensile  Specimens 

Some  previous  investigations  (6)  showed  that  the  anisotropy  in  elastic 
properties  exerts  considerable  influence  on  the  stress  distribution  within 
the  anisotropic  body. 

The  results  of  the  investigation  of  the  influence  of  anisotropic 
mechanical  properties  on  the  stress  distribution  are  presented  here  for  MgO 


10. 


and  A1,0^.  The  principal  reason  for  choosing  these  two  materials  is  the 
availability  of  reliable  elastic  data  in  the  fora  of  compliances.  Also,  the 
crystallographic  symmetries  of  MgO  and  A190^  arc  representative  of  a  large 
number  of  ceramic  materials. 

*». 

The  values  of  the  principal  technv.al  constants  of  MgO  and  A190^  were 
calculated  ^  from  the  compliance  matrices  of  Chung  and  of  Gieske 
respectively. 

Some  selected  values  of  Young's  moduli  of  alumina  are  shown  in  Figs.  4a, 
4b  and  4c.  In  these  polar  diagrams  the  anisotropy  in  the  Young's  moduli  are 
shown  in  the  plane  perpendicular  to  one  of  the  two-fold  axes.  By  a  polar 
rotation  of  approximately  55°  in  alumina,  the  elastic  axes  mav  be  brought 
into  a  position  where  the  difference  between  the  moduli  in  direction  of  the 
hexad  axis  and  one  of  the  diad  axes  is  at  a  maximum.  That  it  is  a  real 
maximum  can  be  seen  from  Fig.  4b  where  the  polar  diagram  of  Young's 
moduli  is  shown  in  the  plane  perpendicular  to  that  of  the  plane  of  Fig.  4a. 
The  other  point  of  significance  may  be  that  where  the  moduli  show  isotropy 
in  the  plane  perpendicular  to  the  tensile  direction.  This  point  is  found  at 
approximately  10°  away  from  the  original  setting  of  axes.  Similarly,  in 
Fig.  4c  the  variation  of  Young's  moduli  is  shown  in  the  plane  perpendicular 
to  the  plane  of  Fig.  4a  and  rotated  by  10°  from  the  original  setting. 

Similar  diagrams  are  presented  in  Figs.  5a  and  5b  where  the  variation 
of  Young's  moduli  of  magnesia  is  shora.  The  maximum  anisotropy  occurs  at 
45°  rotation  between  the  <100  and  <111>  directions.  The  anisotropy  in  the 
plane  perpendicular  to  that  of  Fig.  5a  is  indicated  in  Fig.  5b  where  it  can 
be  seen  that  while  the  absolute  elastic  values  increase,  the  relative 
difference  between  them  reduces.  Based  on  this  set  of  elastic  data,  a  number 
of  representative  orientations  was  chosen  to  keep  the  computational  work 
within  manageable  limits. 
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The  calculated  elastic  properties  and  the  corresponding  orientations 
are  tabulated  in  Table  1. 

It  must  be  nentioned  at  this  stage,  that  the  stress  analysis  of  aniso¬ 
tropic  specimens  of  axisymnetrical  shape  would  require  a  full  three-dinensional 
analysis  to  obtain  meaningful  results.  This,  however,  is  outside  the  capability 
of  the  present  computer  prograa  used  in  these  analyses.  This  work  therefore 
is  restricted  :o  the  analysis  of  flat-hyperbolic  specimen  shapes. 

Figs.  6  ar.d  7  show  the  results  of  analysis  for  the  flat-hyperbolic  speci¬ 
mens  of  alunina  and  magnesium-oxide ,  respectively.  The  stress  distribution  is 
again  shown  as  a  dimensionless  stress  concentration  factor  similar  to  the 
isotropic  case.  It  can  be  seen  that  while  the  general  appearance  of  the 
distribution  curve  is  unchanged,  the  maximum  stress  concentration  ranges  from 
about  1.45  to  approximately  1.32  depending  on  the  selected  combination  of 

elastic  constants. 

Since  the  orientations  were  chosen  to  include  the  limiting  maximum  and 
minimum  elastic  properties,  it  is  probable  that  the  values  of  1.45  and  1.82 
are  the  minimum  and  maximum  concentration  factors  for  all  orientations 
in  alumina. 

Similarly,  the  range  of  stress  concentration  factors  in  magnesium  oxide 
is  found  between  1.45  and  1.65.  It  may  be  significant  that  the  range  of 
stress  concentration  factors  for  alumina  and  magnesia  is  somewhat  in  pro¬ 
portion  to  their  relative  degree  oi  anisotropy. 

E.  Stress  Distribution  in  Anisotropic  Bi-crystal  Tensile  Specimens 

Fc-r  this  part  of  the  analysis  the  bi-crystal  was  assumed  to  be  composed 
by  joining  a  combination  cf  single  crystals  whose  anisotropic  stress  dis¬ 
tribution  was  studied  in  the  previous  section.  This  would  demonstrate  the 
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influence  of  a  boundary  between  two  discretely  adjoined  anisotropic  bodies. 
The  bi-crystal  boundary  teas  chosen  such  that  it  lies  perpendicular  to  the 
tensile  axis  and  is  located  at  the  waist  section. 

In  Table  2  the  combination  of  the  bi-crystals  is  presented  together 
with  a  summary  of  the  "Axial  Notch  Factor"  and  derived  by  FE>1  for  each  of 
the  composing  single  crystals  as  calculated  in  Section  D. 

The  analytical  results  are  presented  in  graphical  ferm  in  Figs.  8  and 
9. 

It  is  easily  observed  that  the  axial  stress  distribution  along  the  focal 
plane  is,  with  the  exception  of  one  case,  unchanged  compared  to  the  axial 
stresses  computed  for  individual  single  crystals.  It  appears  in  general, 
that  the  component  of  stress  perpendicular  to  a  boundary  is  largely  defined 
by  the  elastic  anisotropy  of  the  crystals  located  at  either  side  of  the 
boundary,  and  it  is  influenced  only  to  a  small  extent  by  the  mutual  restraining 
influence  on  each  other  at  the  boundary. 

The  transverse  stress  distribution  in  the  bi-crystals  is  not  as  clear 
as  was  found  for  the  normal  stresses.  For  example,  in  the  bi-crystal 
specimen  A,  the  top  half  of  the  system  showed  a  reduction  in  the  maximum 
value  of  approximately  43%,  allowing  for  a  negative  stress  concentration  at 
the  outer  edges  of  approximately  0.220.  The  lower  half  increased  its 
maximum  by  about  18%. 

A  somewhat  analogous  situation  was  found  in  the  bi-crystal  specimen  D, 
though  the  numerical  values  were  not  so  extreme  as  in  the  case  A  discussed 
above.  In  both  specimens  B  and  C  the  transverse  stress  concentration  was 
found  unchanged. 
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Since  there  were  only  four  combinations  studied  it  could  be  dangerous 
to  attempt  wide  generalizations.  However,  it  appears  that  those  bi-crvstal 
specimens  whose  composing  shear  moduli  are  widely  different  will  display 
changes  in  their  tangential  stress  component,  while  those  bi-crvstals  of 
approximately  the  same  shear  moduli  will  behave  like  two  individual  single 
crystals  with  no  mutual  influence. 

Table  2.  Tabulation  of  FBI  Notch  Factor 


Elastic  Properties 
Identification  from  Table  1 

Notch  Fact 

or 

Upper 

Half 

Lower 

Half 

Upper 

Half 

Lower 

Half 

Isotropic 

Isotropic 

1.45* 

A 

A 

1.71 

B 

B 

1.82 

C 

C 

1.71 

Siugle  Crystal 

D 

D 

1.65 

E 

E 

1.52 

F 

F 

1.47 

G 

G 

1.80 

B 

D 

1.85 

1.61 

B 

G 

1.82 

1.82 

Bicrystal 

B 

F 

1.52 

1.50 

B 

F 

1.79 

1.61 

Same  as  Nueber's  calculated  value  1.454. 


F.  Discussion  of  Results 


A3. 


The  results  show  that  the  elastic  anisotropy  of  the  cater ial  has  an 
effect  on  the  Keuber  notch  factor  for  samples  with  a  hyperbolic  cross  section. 

The  maximum  effect  for  constants  derived  from  alumina,  using  values  for 
maximum  anisotropy,  is  25%  core  than  the  Neuber  value  for  an  isotropic  speci¬ 
men.  The  use  of  two  different  crystals,  top  and  bottom,  does  not  enhance  this 
value,  and  the  two  halves  seem  to  act  independently  of  each  other.  This 
latter  result  indicates  what  seems  to  be  a  deficiency  in  the  model  and  led 
to  the  discovery  of  two  other  results  that  confirm  this  suspicion. 

For  two  cases.  Figs.  9a  and  9d,  the  transverse  forces  do  not  come  to 
zero  at  the  edge  of  the  model.  This  is  a  violation  of  the  known  boundary 
conditions  and  seens  to  indicate  the  presence  of  very  much  higher  stress 
gradients  than  the  present  configuration  of  the  model  can  cope  with.  There 
is  anotuer  indication  that  a  finer  mesh  model  is  needed:  the  integrated 
force,  across  the  focal  plane  for  the  normalized  force  used,  should  equal  one. 
This  condition  is  satisfied  only  for  the  isotropic  case  where  the  FEM  results 
agree  with  the  Neuber  result. 

The  model  used  as  indicated  in  Fig.  2,  is  composed  of  288  triangles  but 
in  the  region  of  the  focal  plane  from  which  the  results  are  obtained  there 
are  only  6  triangles.  It  is  likely  that  this  is  too  coarse  a  mesh  in  this 
region  so  that  the  results  are  affected  by  this  choice  of  model  division. 

The  presence  of  an  enhanced  notch  factor  in  these  results  seems  to 
indicate  that  in  a  refined  model  the  effect  might  be  even  more  pronounced, 
since  it  is  established  that  there  are  stress  gradients  present  that  are 
higher  than  this  model  can  take  into  account. 

G.  Conclusions 

The  finite  element  method  as  applied  to  the  study  of  ceramic  materials 
offers  a  means  of  understanding  the  effect  of  elastic  anisotropy  and  elastic 


16. 


mismatch  on  a  boundary.  The  use  of  the  hyperbolic  shape  adds  a  complicating 
additional  factor  to  the  interpretation,  since  the  shape  effects  the  stress 
concentration  more  than  the  elastic  properties  do. 

In  order  to  evaluate  the  elastic  effects  in  this  configuration  a  finer 
mesh  model  should  be  used. 
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Fig.  2:  The  transverse  half  of  the  F.E.M.  hyperbolic  model. 


Fig,  3a:  Comparison  of  Exact  and  F.E.M.  Isotropic  Solutions  for  Axial 
and  Radial  Stresses  (Flat  Hyperbolic  Specimen). 
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Fig.  3b:  Comparison  of  Exact  ar.d  F.F..M.  Isotropic  Solutions  for  Axial 
and  Radial  Stresses  (Axisymmetr ic  Specimen). 
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Variation  of  Young’s  moduli  for  alumina. 


Fig.  5a:  V  iriation  o£  Young's  moduli  for  magnesia. 
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Fig.  6b: 


•  _c  *.u,  f i -n-_bv»Grbolic  specimens  of  alumina. 
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material*  Identification  B  in  Table  1.  amsotropy) 
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(3-2)  plane 
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Fig .  6d : 


Stress  distribution  of  flat-hyperbolic  specimen  of  single  crystal 
alumina.  Both  top  half  and  bottom  half  are  the  same  anisotropic 
material.  "(3-1)  plane"  refers  to  Identification  E  and  "(3-2) 
plane  refers  to  Identification  F  of  Table  1. 
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Fig.  7:  Stress  distribution  of  the  flat-hyperbolic  specimens  of 
crystal  "(1“2)  plane"  and  "(2-1)  plane"  correspond  to  i 
cations  H  and  J  respectively  in  Table  1. 
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8a:  Stress  distribution  of  the  f ]at-hvperbolic  bicrystal  specimen  of 
alumina.  The  notch  factor  for  the  bicrystal  specimen  exceeds 
slightly  the  values  for  each  material  separately.  The  materials 
are  identifications  B  and  P  of  Table  1. 
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BICRYSTAL  COMBINATION  (B/G) 
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Fig.  8b:  Stress  distribution  of  the  flat-hyperbola  bicrystal  specimen  of 
alumina  using  the  material  properties  of  identifications  B  and  G 
of  Table  1.  The  notch  factor  of  each  separately  is  nearly  the 
same  as  that  for  the  bicrystal. 
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:  Stress  distribution  of  the  flat-hyperbola  bicrystal  specimen  of 
alumina  using  the  material  properties  of  identifications  E  and  F 
of  Table  1,  The  notch  factor  of  each  separately  is  nearly  the 
same  as  that  for  the  bicrystal. 


Preceding  page  blank 


Fig.  8d:  Stress  distribution  of  the  flat-hyperbola  bicrystal  specimen  of 
alumina  using  the  material  properties  of  identifications  B  and 
F  of  Table  1.  The  notch  factor  for  the  bicrystal  lies  within 
the  bounds  for  the  single  crystals.  This  is  in  contrast  to  the 
results  in  Fig.  8a. 
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Fig.  9a:  Transverse  stress  distribution  of  the  flat-hyperbolic  specimens 
of  alumina.  The  materials  arc  identifications  h  and  D  of 
Table  1. 


Transverse  stress  distribution  of  the  flat-hyperbolic  specimens  of 
alumina.  The  materials  are  identifications  R  and  G  of  Table  1. 


Single  Crystal 


Transverse  stress  distribution  o£  the  f lat-hyperbolic  specimens 
of  alumina.  The  materials  are  identifications  I)  and  F  of 
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'  Stress  distributions  in  Al203' 
-Bicrystal  ^Tensile  -Specimens 

J.  E.  GAGORIX,  R.  A.  QUEENEY,  and  H.  A.  McXINSTRY 

D^ostoea  cf  £53025  l5sc!23a  rd  Ceramic  Sotos; 

Tbs  Pts^*os2  Sale  fr.iesly,  Uav«siy  Poi,  Pe=r/tea  IGST2 


Stress  &tis!s&as  ta  A1.0,  bceyslaSs  sta  analyzed  for  cer¬ 
tain  cSseriaaiices  2  cress  tbs  ^eb  toafay  in  sin^ii-sairf 
teasSc  ipecacs  ai!h  the  *rab  teabiy  so.-sul  to  tbs  t  ets-Jc 
axis.  Stresses  wera  found  c sag  tbs  cuatrical  analysis  teeb- 
moys  of  direct  stiflasss  calculations  or  by  fesie  elemfcal  snuf- 
ysis.  5ltn  tbs  ssiseir^  za  elastic  ctsgais  cf  tbs  bccrystals 
stas  considered,  casorisalaSoa  of  priori  pal  axes  yielded 
stress  ccoceairatisns  at  tbs  boundary  as  l^a  as  131  Tbs 
assumption  of  tsnfoaa  stress  suits  ia  bkrystal  racehanical 
response  studies,  therefore,  is  not  generally  acccratt- 


L  Background 

T:t£  role  of  grain  bxabns  ia  detenabscg  tbs  strength 
ptnatiers  of  ptfycrygaHiag  sands  has  bsca  ismts- 
gated  extensively.  For  polycrysmlline  caffes  of  tow  or  vatstb- 
iag  ductility,  fracture  strcrgtbs  are  affected  strong!..  k  tbe 
presence  of  grab  banduics.  Many  of  tbs  known  fads  abort 
tbs  effects  of  grab  boiadaries  ca  tbs  mechanical  properties  of 
ceramics  are  reviewed  b  Ref.  1. 

Oae  method  of  studying  the  strength  characteristics  of  grain 
boundaries  is  through  the  use  of  k crystal  specimens.  Ia  prin¬ 
ciple,  the  rnisorienlalsen  of  the  adjoining  crystals  and  the 
orientation  o!  the  boundary  can  be  coat  rotted  to  provide  a  ra¬ 
tional  program  of  strength  testing.  The  ctfiniuoa  of  bicrystal 
geometry*  and  characterization  of  fracture  strength  of  MgO 
bicrystals  was  studied  by  Lange  and  Btics-em  ~*  Similar 
studies  have  been  conducted  on  NaCi  bicrystals.*  The  high- 
temperature  creep  resistance  of  grain  tour.carics  in  pressure- 
sintered  A!.0j  bicrystals  has  been  investigated.* 

A  common  feature  of  these  studies  was  the  assumption  of 
uniform  stress  states  at  and  near  the  bicrystal  boundary.  That 
this  assumption  is  not  exact  can  be  seen  qualitatively  by  con¬ 
sidering  a  soft  material  joined  :«  a  stiffer  one  with  a  tensile 
load  applied  normal  to  the  boundary.  The  larger  tensile 
strains  in  the  softer  material  iv:il  produce  larger  transverse 
strains  as  well,  assuming  roughly  equal  Poisson’s  ratio  effects. 
A  transverse  stress  field  near  the  boundary,  therefore,  is 
necessary  to  preserve  compatibility  cf  deformation.  This  re¬ 
sult  is  not  predicted  when  3  continuously  uniform  stress  field 
at  the  boundary  is  assumed. 

Bicrystal  specimens  can  be  treated  similarly  because  the 
misorientation  of  the  components  generally  results  in  a  change 
in  clastic  properties  along  a  giver,  direction  at  the  boundary. 
The  non-nominal  stress  slates  developed  during  loading  should 
be  defined  before  mechanical  strength  measures  can  be  speci¬ 
fied  accurately. 


II.  Method  of  Analysis 

A  closed-form  analytical  solution  to  the  stress  field  near  the 
transverse  grain  boundary  in  a  bicrysta!  tensile  specimen  is  a 
formidable  problem  because  of  the  finite  width  of  real  speci¬ 
mens.  Expanding  the  specimen  geometry  to  two  joined  half¬ 
spaces  might  simplify  the  analytical  solution,  but  this  approach 
would  not  be  applicable  to  the  finuc-uidtli  real  geometry. 
Furthermore,  the  most  extreme  values  of  the  'tress  distribu¬ 
tion  may  be  developed  at  the  specimen  edges.  This  informa- 
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isca  would  be  of  prime  ixpxteete  ia  atearical  strength 
studies. 

To  efreumveut  these  problems,  the  numeric?!  calculation 
rmhai  of  fiaite^iement  analysis  was  used.  Toe  continuous 
structure  (bicrysta!)  is  subdivided  into  a  set  of  clastic  fiurie 
elements  (triangles  in  the  present  work).  The  triangles  are 
joined  discreetly  at  the  nodes.  A  complete  stiffness  matrix 
for  the  structure  can  be  assembled  by  adding  the  stiffnesses 
of  each  element  at  the  nodes.  The  stiffness  matrix  relates 
forces  on  the  nodes  to  nodal  displacements  by: 

[F]=[K][«]  CD 

By  inverting  the  stiffness  matrix,  the  nodal  displacements  can 
be  related  to  the  nodal  forces  by: 

[a]=[K]![F]  (2) 

The  stress  in  each  clement  can  be  found  from  the  nodal  dis¬ 

placements  according  to: 

W=[S][S1  (3) 

Thus,  the  components  of  the  stress  tensor  [a]  is  any  element 

can  be  found  from  the  stiffness  matrix  for  the  element,  [SJ, 
and  the  nodal  displacements  calculated  in  Eq.  (2).  An  am¬ 
plification  of  the  finite-element  method  is  given  by  Zicnkicv.icz.* 
The  method  of  analysis  has  bcc-n  applied  to  the  intcrfaciai 
stress  states  in  bicrystals/ 

III.  Bicrystal  Geometry’  and  Properties 

The  specification  of  crystalline  misorientation  across  the 
grain  boundary  is  shown  in  Fig.  1,  in  which  the  unprimed  axes 
are  the  principal  axes  of  the  unit  cell:  coincident  with  these 
are  the  underlined  Euler  axes.  The  primed  axes  define  rota¬ 
tions  from  the  principal  frame  through  Euler  angles  a  and  p, 
as  indicated. 


Received  April  5,  1971;  revised  copv  received  June  1 4, 
1971. 

Supported  in  part  by  the  United  States  Office  of  Naval  Re¬ 
search  under  Grant  No.  N0014-G7-A-03S5-0004. 


64. 

Journal  of  The  American  Ceramic  Society — Gagorik  el  aL  Vol.  54,  No.  12 


FIp.  —  Tcn-TIc  »pec- 

inxn  liinfi"  detrapo- 
iiiisn  into  finile  clc- 
cacat*;  main  toBsiaij 
located  at  dated  line. 


VK!  9  »7 

rtor  the  purposes  off  the  peasant  sandy.  a  pfene-saass, 
straight-saded  tsssSe  specimen  OiHJD  ta.  tag:  C  450  Ea  wife. 
sa3  0.033  a  eras  casta  Tra  spjdaa  aaS  its  tek- 
dhma  into  tri jeguta-  ttetsa  are  shown  in  Fsg  2  Thus  span* 
men  aryl  fibs  boundary  candhaccs  imposed  on  it  cerrespced  to 
a  p3aas  stress  ccnf^ururion.  Tbe  eexe  of  a  ITIkSaes 
r-^-A  «as  arrived  ai  by  trial  and  cnar.  Tee  goal  aas  to 
ariaeve  a  spedmea  ia  which  a  nominal  tzslmxa  Cress  state 
was  rcalired  near  tbe  sperimta  ends:  Accsa&ad  cataihs- 
tica  error,  maaitorad  ia  the  form  cf  departure  fraa  load 
c^a&ssa.  was  tnamlafaed  at  <1*J. 

Tbe  etoea  stiffnesses  [S]  were  forcisri  fears  fte  sig;*- 
crystal  dasric-censaam  stiffness  compe-eds  C  !?  C-»>  ate 
the  principal  axes  compcssas  c:  tbc  stifK-ss  cocSdcats,  the 
componcsss  ia  a  new,  primed  tads  frame  ate: 

(4) 

with  tbe  summation  convention  observed.  Tbe  transformation 
coeSudcnts  o„  are  the  direction  cosines  between  the  new  and 
old  axes  frames.  The  tschcrmsl  stiff  nesses  used  for  ALOi  were 
those  of  Aleksandrov  aad  Ryabova.* 


IV.  Results  and  Discussion 

Changes  ia  stress  states  caused  by  ralsori entations  between 
the  upper  and  lower  halves  of  the  bicr-stal  were  examined. 
To  dismiss  these  changes  over  the  nominal  state,  the  stress 
concentration  factor,  X„  may  be  used,  where 

.  actcal  axial  stress  state 
=  nominal  axial  stress  state 

Since  the  total  principal  stress  state  determines  mechanical 
response  in  general,  the  development,  with  misorientation.  of 
transverse  stresses  should  also  be  characterized.  Since  there 
is  co  nominal  or  average  transverse  stress,  the  usual  definition 
of  stress  concentration  factor  cannot  be  used.  Instead,  a 
transverse  stress  concentration  factor  is  defined  as  the  ratio 
of  Uie  transverse  stress  to  the  nominal  axial  stress. 

Results  obtained  are  shown  in  Figs.  3  and  4.  Two  types  of 
misorientations  were  studied. 

Case  A:  The  plane  of  the  specimen  has  its  normal  in  the 
3-direction  of  the  unit  cell  with  tension  applied  in  the  1-direc¬ 
tion.  The  bottom  half  of  the  bicry-tal  remains  stationary  for 
all  calculations,  whereas  the  top  half  is  rotated  through  various 
a  angles,  with  the  angle  3  held  constant  at  zero.  Figure  3(A) 
indicates  that,  for  a  40°  a  rotation  of  the  upper  half,  an  axial 
stress  concentration  factor  of  stl.4  is  ach.evcd  at  the  specimen 
side  edge.  This  figure  also  indicates  a  monotonically  increas¬ 
ing  stress  concentration  factor  as  misorientation  is  increased. 
Figure  3(B)  demonstrates  the  tendency  of  the  stress  concen¬ 
tration  factor  to  return  to  unity  as  one  moves  away  from  the 
grain  boundary,  as  expected.  Figure  3(C)  demonstrates  the 


ggvdopcnent  d  transverse  cresses  with  increasing  nsscrieaa- 
tica.  Toe  transverse  stresses  decease  n.a£y  with  £s&ace 
from  tbe  grim  boundary  and  do  not  vary  appreciably  abb 
dfstnnrie  a tag  tbe  interface. 

Case  3:  Tbe  top  half  of  tbe  specimen  was  originally 
orient td  as  in  Case  A.  Tbe  bcaon  half  is  oriented  wfth  3- 
dlrettiaa  ia  tbe  tensSe  dacha,  I-direerien  in  tbe  transverse 
sense,  and  2-dbectka  normal  to  tbe  piane  off  tbe  specimen. 
Tbe  bottom  half  was  hric  fixed  Its  all  calculations,  wbireas 
tbe  epper  fealf  was  rotated  as  instated  in  Fig.  4. 

Aadal  stress  ccoreairaiioa  faeries  as  large  as  13  were 
found  at  tbc  specimen  edge  (rig.  4(A)).  Figure  4(3)  EnS- 
cates  that  these  valaes  tend  mriddy  toward  eaity  away  from 
tbe  grain  bxnffsy-  Large  transverse  stresses  wrath  were 
igjiy  tocalrred  at  tbc  grain  boundary  also  devdoped  (Fig. 
4(0);  they  «5d  not  change  appreciably  acres  tbe  specimen 
width. 

Tbtse  resells  indstale  that,  as  a  generality,  raisorieaEEitsa 
between  tbc  halves  of  an  ALO,  bkrysta!  tends  to  develep  am 
severe  stress  states  at  tic  grain-bssndary-spechsen-edg2  ia- 
terseclraa  than  were  previses?,  tbrcgbl  to  exist.  Tic  presr=t 
calculations  found.  for  certain  selected  rcsoricataiams, 
stresses  50%  iagber  than  those  predicted  by  tbe  assumption 
cf  average  comma!  axel  stress  for  tbe  specimen.  Transverse 
stresses,  winch  cataraHy  affect  she  saedmmtai  response  off 
tbe  fcaystal.  were  also  sound  to  be  present  and  large  enough 
to  be  sometimes  non-negligicic. 

Flares  3  and  4  indicate  that  tbe  calculated  axial  stresses 
are  not  sdenucal  at  the  interface,  ia  apparent  violation  of  force 
cfaBnan  Tins  discrepancy  is  an  intrinsic  characteristic 
cf  the  numerical  techaanse  used.  Tbe  sinite-ctancut  rrecccd 
averages  the  stresses  over  tbe  entire  dement  and  associates 
tbe  average  value  with  tbe  dement  centroid.  Although  tbe 
particular  programing  arrangement  used  in  the  present  study 
ensured  nodal  equilibrium  to  within  1.0%.  no  such  demands 
can  be  imposed  on  these  average  values  for  tbe  entire  dement. 
Tbe  use  of  element  average  stresses  coses  no  problem  unless 
the  analysis  is  concerned  with  abrupt  changes  in  geometry  or 
elastic  properties.  A  finer  mesh  choice  at  the  interface  would 
undoubtedly  indicate  equilibrated  axial  stresses  there.  How¬ 
ever.  the  prcse.it  study  was  designed  to  draw  attention  to  the 
necessity  of  more  accurate  stress  analyses  for  anisotropic 
bodies,  not  to  characterize  given  experimental  specimens 
thoroughly. 

it  is  cot  possible  to  draw  firm  conclusions  about  the  relative 
effects  of  variations  in  specific  clastic  constants  at  the  interface 
on  the  stress  distribution  there.  Voile  variations  in  eg.  the 
axial  stiffness  differences  across  the  boundary  are  bring  in¬ 
vestigated.  the  other  elastic  constants  are  also  changing,  and 
the  stress  distribution  (Eq.  (3))  depends  on  all  of  these. 

In  view  of  the  results  of  the  present  study,  experiments  in¬ 
volving  the  stress-dependent  properties  of  fcicryslal  grain 
boundaries  must  include  an  analysis  of  the  stress  state  that  ac¬ 
counts  for  the  rapid  change  in  elastic  properties  at  the  inter¬ 
face.  This  requirement  would  apply  to  all  symmetry  classes 
of  crystal  structure,  not  the  trigonal  class  tc  which  AI.O,  be¬ 
longs  only.  The  stress  analysis  would  be  viewed  in  the  same 
characterization  context  as  determination  of  the  crystallo¬ 
graphic  orientations  cf  the  crystal  halves. 

V.  Summary 

The  abrupt  change  in  clastic  properties  at  the  grain-bound¬ 
ary  interface  in  bicrystal  tensile  specimens  (\iO,  was  in¬ 
vestigated)  prodaced  changes  in  the  nominal  stress  state  there. 
These  new  stresses  are  highest  at  the  specimen  edge,  localized 
near  the  boandary,  and  vary  with  crystallographic  mismatch. 
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ABSTRACT 

The  elastic  behavior  of  a  ceramic  body  has  been  modeled  by  use  of 
the  finite  element  method.  The  model  simulates  a  set  of  nineteen  hexagonal 
crystals  whose  elastic  and  thermal  properties  can  be  independently  varied. 
Several  different  configurations  of  the  elastic  properties  of  the  crystals 
were  evaluated  under  conditions  of  isostatic,  uniaxial  and  thermal  loadings. 
The  effect  of  anisotropy  shifts  the  position  where  the  maximum  shear  acts 
in  the  model,  under  isostatic  and  thermal  loadings. 

INTRODUCTION 

Tha  finite  element  method  of  analysis  has  been  groxjing  rapidly  in 
application  to  widely  varying  subjects:  from  earth  dams,  to  whole  steam¬ 
ships  to  pressure  valves  (Refs.  1,2).  The  present  study  may  be  an  applica¬ 
tion  to  the  smallest  objects  yet  studied:  the  grains  of  a  ceramic  body. 

The  model  for  study  was  taken  to  be  the  simplest  two-dimensional 
representation  of  a  ceramic  body.  Buessem  and  Lange  (Ref.  3)  proposed  a 
hexagonal  model  in  which  the  forces  due  to  thermal  expansion  anisotropy 
might  be  considered.  In  his  discussion  of  the  microstructure  of  ceramics 
p.  409  (Ref.  4),  Kingery  states,  "if  we  restrict  the  structure  to  one  with 
the  simplest  corners  with  three  grains  meeting,  the  average  polygon  must 
be  a  hexagon."  The  model  designed  for  this  study  consists  of  nineteen  hexa¬ 
gon-shaped  crystals  as  shown  in  Fig.  1.  Each  of  these  hexagons  can  be  given 
specific  anisotropic  elastic  constants  and  anisotropic  thermal  expansion 
coefficients.  In  this  way  the  model  can  simulate  the  random  orientation  of 
crystals  in  a  polycrystalline  ceramic  body. 

This  model  was  further  subdivided  into  366  triangles  with  199  nodes  as 
shown  in  Fig.  2,  in  order  to  be  compatible  with  the  computer  program  available. 

The  computer  program  was  written  by  Dahl  (Ref.  5)  and  has  been  modified 
by  Gagorik  (Ref.  6)  and  T.  Kovacs  (Ref.  7). 


Alumina  Model 


The  elastic  constants  of  alumina  were  used  to  obtain  the  values 
assigned  to  each  of  the  nineteen  hexagonal  crystals  of  the  model.  The 
elastic  compliances  used  were  those  measured  by  J.  Geiske  (Ref.  8).  The 
values  used  were  converted  into  the  English  system  of  units  as  are  given 
in  Table  I. 


Table  I 

2 


Elastic  Compliances 
Si  i  =  0.01620  S24 

512  =  0.00480  S22 

5 1 3  =  0.00264  S5s 

514  =  0.00327  Sss 

S33  »  0.01500  S66 

Si --  0.00478 


(inches)  / pound 
=  -Sii, 

-  Si  1 
=  S44 

=  -2S 14 

=  2  (Sn  -  Si  2) 


These  elastic  compliances  were  first  transformed  by  a  computer  program 
written  by  J.  Gagorik  (Ref.  6)  based  on  the  method  of  Lieberman  and 
Zirinsky  (Ref.  9)  which  performs  the  equivalent  of  a  rotation  of  the  crystal - 
Then,  these  compliances  for  the  rotated  crystal  were  converted  into  technical 
constants.  For  an  anisotropic  crystal,  technical  constants  are  defined  by 
Lekhnitskii  (Ref.  10)  to  be  consonant  with  the  usual  elastic  constants; 
Young's  modulus,  E,  the  shear  modulus  G,  and  Poisson's  ratio,  V. 

For  the  two-dimensional  plane-stress  problem  considered,  the  technical 
constants  are  obtained  from  the  value  of  the  rotated  elastic  compliances  by 
the  following  relationships: 
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E>;  =  1/Sn 
Ey  =  I/S22 

\y  =  -Sl2/Sll 

v  =  -S12/S22 

G  =  1/S66 
xy 


Hooke's  law  in  plane  stress  is  then  expressed  as 
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Configurations 

A  set  of  six  different  configurations  was  processed  using  four 
different  types  of  leading. 

In  the  first  configuration,  each  of  the  nineteen  crystals  was  given 
isotropic  constants  corresponding  to  the  orientation  of  each  crystal  with 
the  c-axis  perpendicular  to  the  plane  of  the  model. 

The  second,  third  and  fourth  configurations  were  obtained  by  assigning 
to  the  center  crystal  technical  constants  corresponding  to  the  rotation  of 
an  alumina  crystal  about  an  axis  parallel  to  the  x-axis  by  amounts  of  40, 

90  and  140  degrees  respectively.  The  rest  of  the  constants  were  given  the 
same  isotropic  values  as  in  the  first  configuration. 

The  fifth  configuration  was  composed  of  three  different  sets  of  con¬ 
stants.  The  center  crystal  was  oriented  so  that  the  c-axis  lies  in  the  plane 
of  the  paper  and  points  in  the  x-direction.  The  a-axis  was  supposed  to 
have  been  placed  in  the  y-direction.  However,  a  value  for  the  direction 


corresponding  to  the  maximum  anisotropy  was  selected  instead.  The  regain¬ 
ing  eighteen  crystals  were  divided  in  such  a  way  that  for  the  top  half  of 
the  crystals  all  had  the  same  orientation,  and  the  bottom  half  crystals 
all  were  given  constants  corresponding  to  a  ninety  degree  rotation  from  the 
other  half  of  the  crystals.  The  constants  for  these  two  halves  were  selected 
to  correspond  to  the  orientation  found  by  Kovacs  (Ref.  7)  to  possess  the 
maximum  anisotropy  for  alumina.  This  corresponds  to  a  rotation  of  the 
c-axis  of  the  crystal  about  an  axis  perpendicular  to  an  a-axis.  The  rotation 
for  maximum  anisotropy  was  fifty-five  degrees.  In  the  model,  constants  were 
used  corresponding  to  a  rotation  of  fifty  degrees. 

The  sixth  configuration  was  an  arrangement  designed  to  represent,  as 
well  as  possible,  a  polycrystalline  ceramic  body  with  anisotropic  constituents. 
The  constants  were  chosen  to  produce  a  random  arrangement  of  the  crystals. 

The  center  crystal  was  fixed  in  its  orientation  so  that  the  c-axis  was 
placed  in  the  x-direction  and  the  a-axis  in  the  y-direction.  A  two-angle 
rotational  matrix  with  ten-degree  intervals  in  both  angles  was  used  as  the 
population  from  which  to  choose  the  appropriate  constants  randomly. 

Loadings 

Each  of  the  six  configurations  was  subjected  to  four  different  loadings. 
The  first  loading  corresponded  to  an  isostatic  loading.  Force  components 
were  placed  on  each  external  node  in  such  a  way  that  the  resultants  would 
point  to  the  center  of  the  model.  This  type  of  loading  produces  some  shear 
in  the  outermost  elements  but  is  evened  out  before  reaching  the  second 
layer  of  crystals. 

The  second  and  third  loadings  were  uniaxial  tensions,  applied  at  the 
external  nodes.  In  the  second  type  the  forces  act  vertically  in  the  y~ 
direction  and  in  the  third  type,  the  forces  act  horizontally  in  the  x- 
direction.  Hence,  the  y-direction  loading  is  uniform.  The  x-direction 
loading  is  slightly  non-uniform.  In  the  fourth  type  of  loading,  the  forces 
generated  by  the  thermal  expansion  of  the  crystals,  with  each  expansion 
determined  by  its  appropriate  anisotropic  thermal  expansion  coefficients, 
vere  placed  on  the  external  nodes.  The  system  was  then  allowed  to  adjust 
internally  under  these  forces.  For  the  isotropic  case,  it  was  found  that 
the  forces  were  not  identical  with  those  used  for  the  isostatic  case.  The 
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forces  were  calculated,  however ,  so  that  no  shear  strains  were  generated 
at  any  point  within  the  nodel. 

Finite  Element  Method 

The  finite  element  method  provides  a  means  of  obtaining  the  solution 
of  elastic  problems  which  are  not  directly  soluble  by  standard  elastic 
mathematical  analysis.  The  structure  to  be  analyzed  in  divided  into  a  set 
of  small  units,  usually  triangles,  to  which  the  standard  elastic  theory 
can  be  applied.  Each  triangle  is  linked  with  the  triangles  surrounding  it. 
The  displacements  of  the  connecting  nodes  then  must  be  related  to  the  other 
in  a  compatible  way,  and  the  computer  program  finds  the  answer  as  a  solution 
of  simultaneous  equations. 

The  time  required  for  computer  solution  varies  with  the  parameters 
assigned  to  the  model  and  ranges  from  60  seconds  to  120  seconds.  The  lower 
time  is  the  more  common.  The  program  uses  a  Gauss-Seidel  iterative  method 
for  solution,  and  there  is  a  variation  in  the  number  of  iterations  required 
to  reduce  the  error  to  an  acceptable  level. 

Representation 

The  figures  (Figs.  3-10)  which  present  the  results  of  these  analyses 
show  on  one  page  all  six  configurations.  For  each  different  type  of  loading 
there  are  two  pages;  one  showing  the  principal  stresses,  and  the  other  show¬ 
ing  the  maximum  shear  stress. 

In  the  figures  presenting  the  principal-stress  data,  the  length  of  the 

half  arrow  represents  the  magnitude  of  the  maximum  tensile  stress  at  the 

center  of  each  element,  and  the  magnitude  of  the  second  principal-stress  is 

indicated  by  a  line  perpendicular  to  and  crossing  it.  The  direction  in 

which  these  stresses  act  is  also  indicated.  In  the  figures  showing  the 

maximum  shear  stress  (t  ),  the  direction  indicated  is  fortv-five  degrees 

max  • 

to  the  principal  stresses,  and  the  magnitudes  of  both  members  of  the  cross 
are  equal.  In  order  to  show  the  variations  in  T  more  clearly  for  the 
isostatic  and  thermal  loading  cases,  a  scale  factor  was  used  to  increase 
the  size  by  about  16%  over  that  used  for  the  principal  stresses. 


The  computer  output  consists  of  a  tabulation  of  the  stress  and  strain 
components  at  the  center  of  each  triangle. 

DISCUSSION 


Isostatic  Loading 

One  of  the  interesting  features  is  evident  in  comparing  the  t 
plots.  Pigs.  4b.  c  and  d  for  the  40°,  90°  and  140°  rotation.  The  maximum 
shear  stress  lies  outside  the  central  crystal  for  rotations  of  40  and  90 
degrees.  It  lies  inside  for  the  140°  rotation.  The  effect  of  the  aniso¬ 
tropy  is  to  vary  the  position  in  the  sample  where  "nax  will  have  its 
largest  value. 

In  the  bicrystal  with  a  central  inclusion.  Fig.  4e,  the  distribution  of 

T  has  a  semblance  of  a  2-fold  symmetry  axis  outside  the  central  crystal, 
max 

The  greatest  shear  lies  inside.  Its  magnitude  is  approximately  10%  of  the 
principal  stresses.  The  central  crystal  is  anisotropic,  and  this  adds  an 
additional  complication  to  the  interpretation  of  the  source  of  the  shear 
intensif icati on . 

The  random-orientation  model  gives  a  maximum  shear  stress  near  a  corner 
of  the  central  crystal.  The  magnitude  of  the  Tmax  is  a  little  more  than  6% 
of  the  principal  stresses. 

Uniaxial  Loadings 

Two  different,  compressional  uniaxial  loadings  were  used,  since  their 
action  at  different  faces  of  a  hexagon  would  be  different.  The  figures 
(Figs.  5, 6,7,8)  offer  little  help  in  an  interpretation  of  the  effects  of 
anisotropy.  The  lines  of  principal  stress  align  themselves  with  the  external 
loads.  There  is  some  evidence  of  a  slight  turning  but,  as  read  from  the 
computer  output,  the  variation  of  the  angle  of  principal  stress  varies  by 
not  more  than  ±  2  degrees.  (It  is  interesting  to  note  that  the  direction 
of  principal  stress  coincides  with  the  principal  strain  in  the  isotropic 
media  and  deviates  a  little  in  anisotropic  media).  By  extrapolating  the 
principal  stresses  to  the  faces  of  the  central  crystal,  a  maximum  variation 
of  6%  in  the  different  models  is  observed  in  the  face  normal  stress  and 
the  face  shear  stress. 


The  changes  in  the  normal  force  on  faces  one  and  four  for  the  case 
of  vertical  uniaxial  loading  can  be  obtained  from  the  computer  output* 

The  stresses  normalized  in  terms  of  the  applied  stress  for  the  40%  90° 
and  140°  rotations  are  1.5,  0.3  and  -5.0  percent  respectively.  The  effect 
of  the  anisotropy  is  sufficient  to  change  the  sign  of  the  stress  from  tension 
to  compression. 

In  an  attempt  to  gain  more  information  about  the  uniaxial  loadings,  a 
different  method  of  presenting  the  data  has  been  shown.  The  results  shown 
in  Figs.  11-19  xri.ll  be  discussed  under  the  heading  of  SYMAP  contour  maps. 

Thermal  Loading 

For  the  isotropic  case,  thermal  loading  of  the  model  produced  absolutely 
no  internal  stresses,  (Fig.  9a).  With  the  rotation  of  the  central  crystal 
about  an  axis  parallel  to  the  x-axis,  interesting  differences  can  be  noted. 

The  principal  stresses  maximize  inside  the  center  crystal  for  the  40°  rota¬ 
tion  case  and,  as  can  be  seen  in  Fig.  9b,  there  is  a  tendency  for  uniaxial 
stressing  on  the  left,  outside  the  center  crystal  boundary.  They  maximize 
outside  and  to  the  left  of  the  center  crystal  for  the  90°  rotation,  (Fig.  9c). 
For  the  140°  rotation,  (Fig.  9d),  the  principal  stresses  become  almost  uni¬ 
axial  on  the  right  side  of  the  crystal  just  inside  the  boundary.  The  tra¬ 
jectories,  as  indicated  by  the  arroxjs,  shoxj  that  in  addition  to  the  anisotropy, 
the  geometrical  shape  of  the  boundary  also  has  a  strong  effect.  The  high 
T  near  the  boundaries  of  the  center  crystal  emphasizes  the  alteration  of 
the  distribution  of  stresses  as  the  anisotropic  character  of  the  crystal  is 
changed,  (Figs.  10b, 10c  and  lOd).  The  distributions  are  clearly  different 
in  the  three  cases. 

The  normal  and  tangential  stresses  on  each  face  of  the  center  crystal 
were  calculated  by  extrapolating  the  stresses  to  the  boundary.  Parabolic 
extrapolation  functi.ons  were  used  to  find  the  best  fit  from  four  nearest 
layers  both  inside  and  outside  the  center  crystal.  The  inner  and  outer 
extrapolations  do  not  always  agree.  For  Table  II,  which  gives  the  values 
for  these  three  rotation  models,  the  maximum  deviation  for  inner  and  outer 
extrapolations  was  approximately  50%  with  the  average  deviation  about  10%. 

For  the  isotropic  case  these  stresses  xrould  be  zero. 
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Table  II 

The  normal  and  tangential  stresses  (psi)  on 
the  faces  of  the  center  crystal 


40  degrees  90  degrees  140  degrees 


Face 

fr 

He 

^T 

1 

3800 

-70 

4400 

-65 

-1466 

0.2 

2 

15015 

+2810 

14876 

2900 

-2670 

-994 

3 

10430 

-1475 

10350 

-1681 

-6200 

1800 

4 

9010 

+35 

9700 

24 

4074 

-70 

5 

10400 

1500 

10360 

1720 

-6150 

-1790 

6 

14940 

-2794 

14810 

-2905 

2650 

1000 

The  values  for  tha  40°  and  90°  rotations  are  much  more  similar  than 
the  diagrams  indicate.  The  140°  rotation  is,  however,  markedly  different. 
Whereas  the  first  two  leave  the  center  crystal  in  tension,  the  third  case 
has  it  in  a  somartiat  more  complicated  state,  with  four  faces  experiencing 
compression  ..nd  Lwo  faces  under  tension.  As  indicated,  these  stresses  should 
not  be  regarded  as  having  precision.  Since  some  of  the  inner  and  outer 
extrapolation  values  deviate  by  as  much  as  50%,  it  must  be  assumed  that  the 
model  is  not  sufficiently  .*  abdivided  to  take  into  account  the  rapidly 
changing  stress.  With  the  stress  gradient  so  high  in  the  region  of  interest, 
a  new  model  would  be  indicated  as  necessary.  Unfortunately,  time  and  funds 
do  not  permit  such  an  extension. 

Although  the  deviations  of  the  stresses  for  the  maximum  anisotropy  case 
vary  on  an  average  of  45%,  all  of  them  have  the  same  sign.  This  is  not  the 
situation  for  the  random  orientation  model.  Here,  the  deviations  may  be  as 
large  as  500%  and  the  extrapolated  values  are  of  opposite  sign.  Clearly,  a 
finer  mesh  model  is  needed  for  this  case.  The  normal  and  tangential  stresses 
range  from  plus  or  minus  zero  to  thirty  thousand  psi.  It  is  unfortunate 
that  our  present  model  does  not  give  reasonable  values  for  this  nearest 
approach  to  the  real  ceramic  body. 

SYMAP  Contour  Hap 

Another  method  for  presenting  the  data  from  the  FEM  program  is  to  make 
use  of  an  available  computer  program  SYMAP  (Ref.  11)  which  prints  out  an 
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interpolated  contour  map. 

The  computer  over-prints  characters  to  vary  the  darkness  of  the  print¬ 
ing,  and  thus  a  contour  map  is  developed  for  visual  inspection.  The  results 
for  the  uniaxial  loadings  of  the  rotated  center  crystal  are  presented  in 
Figs.  11-16.  The  maximum  values  do  not  appear  in  these  figures,  because  only 
the  region  of  and  near  (four  layers  out)  the  center  crystal  are  presented. 
Where  the  maxima  lie  within  this  region,  the  figures  are  easily  interpreted 
since  they  are  the  plot  of  a  single  variable  only.  However,  as  an  aid  to 
interpretation,  contour  lines  separating  the  different  stress  fields  have 
been  hand-drawn. 

The  data  plotted  in  Figs.  11-16  are  the  maximum  principal  stress,  <Ji. 
Since  the  uniaxial  loading  was  compressions!,  the  second  principal  stress  will 
be  negative  and  of  larger  magnitude.  Thus  the  Oj  plotted  here  will  be  a 
sensitive  measure  of  the  deviation  from  a  uniform  uniaxial  stress  distribution. 
For  a  completely  uniaxial  distribution,  Oj  should  be  zero. 

A  scale  that  can  be  used  to  describe  the  maps  was  obtained  by  dividing 
the  calculated  value  for  0*  by  the  absolute  value  of  a 2.  This  will  produce 
a  ratio  of  values  from  -0.22  to  +0.08  on  a  scale  of  8  as  listed  in  Taole  III. 

In  Figs.  11-16,  the  scale  number  is  produced  at  the  center  of  each  triangle 
in  the  finite  element  model. 


Table  III 

Scale  for  SYMAP  Contour  Plots 


Scale 

Number 

L 

1 

2 

3 

4 

5 

6 

7 

Cfl 

-0.28 

-0.22 

-0.18 

-0.12 

-0.08 

-0.02 

+0.02 

+0.08 

M 

Figs.  11,  12  and  13  correspond  to  the  plots  in  Figs.  5b, c  and  d,  for 
the  rotation  of  the  center  crystal  by  40°,  90°  and  140°  respectively. 
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For  the  40°  rotation.  Ci  has  its  largest  magnitude,  -0.28,  to  the 
right  and  to  the  left  of  the  center  crystal  in  two  places  adjacent  to  the 
corners.  The  pattern  of  the  stress  distribution  for  the  90°  rotation  is 
similar.  The  appearance  of  the  140°  rotation  is  considerably  different. 

The  maximum  compression  is  shifted  to  the  top  and  bottom  of  the  center 
crystal.  The  compression  at  the  center  of  the  center  crystal  changes  from 
-0.12  to  -0.18  to  -0.22,  respectively,  for  the  three  rotations. 

The  uniaxial  vertical  loadings  show  a  considerably  different  pattern. 
Figs.  14,  15  and  16  correspond  to  the  plots  in  Figs.  7b, c,  and  d  for  the 
rotation  of  the  center  crystal  by  40°,  90°  and  140°,  respectively.  The  O i 
is  slightly  in  tension,  +0.02  for  the  40°  and  90°  cases,  but  it  is  com- 
pressional  for  the  140°  case,  -0.08.  There  is  a  corresponding  build-up  of 
compressional  stress,  top  and  bottom  in  these  cases,  perpendicular  to  the 
uniaxial  stress,  but  it  achieves  a  maximum  value  of  only  -0.12  for  the  40° 
case. 

Off  Center  Models 

A  slightly  different  set  of  models  is  presented  in  Fig.  17.  The  crystal 
to  the  right  of  the  center  one  has  been  rotated  40°,  90°  and  140°,  re¬ 
spectively,  in  Figs.  17a, b  and  c,  and  subjected  to  isostatic  loading.  Both 

T  and  the  principal  stresses  are  presented  on  the  same  page.  The  T 
msx  max 

plots  show  much  more  variation  than  do  the  plots  of  the  principal  stresses. 

The  distribution  of  the  maximum  shear  stresses  indicates  a  rather  profound 
effect  of  anisotropy.  The  larger  values  of  the  shear  for  the  40°  case  are 
outside  of  and  to  the  left  of  the  rotated  crystal,  with  a  considerable  build¬ 
up  on  the  slant  faces. 

The  90°  case  shows  a  similar  distribution,  though  the  magnitude  of  the 

T  is  not  as  large, 
max  ° 

The  distribution  of  T  in  the  140°  case  is  markedly  different,  with 

max 

the  larger  values  lying  inside  the  boundaries  of  the  rotated  crystal.  Similar 
behavior  was  observed  in  Figs.  4b, c  and  d,  when  the  rotated  crystal  was  the 
center  crystal. 


i 


it 
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Simplified  Models 

The  results  of  the  finite  element  analysis,  are  plotted  in  a  slightly 

different  form  in  Figs.  18c,  d  and  19.  Each  principal  stress  and  T 

max 

are  plotted  separately.  The  first  principal  stress  (maximum  tension)  is 
plotted  on  the  left  as  a  single  half  arrow  at  the  center  of  each  triangle 
of  the  model.  The  length  of  the  arrow  is  proportional  to  the  magnitude  of 
the  stress  it  represents.  The  magnitude  of  the  maximum  arrow  length  is 
indicated  in  this  figure  for  each  stress.  The  second  principal  stress 
(maximum  compression)  is  plotted  similarly  in  the  center,  and  the  maximum 
shear  stress,  T  x>  is  plotted  on  the  right.  Only  isostatic  loading  is 
considered. 

In  the  case  of  Model  A,  a  hard  isotropic  inclusion  is  inserted  into 
the  system.  The  Young's  modulus  of  the  center  crystal  is  ten  times  larger 
than  the  surrounding  crystals.  The  stresses  on  all  faces  are  essentially 
equal  and  raised  by  the  factor  1.36  over  that  of  the  completely  homogeneous 
model.  At  the  boundary  the  stress  is  seen  to  be  predominantly  radial.  The 
extrapolated  normal  stress  at  the  boundary  is  consistently  3%  higher  inside 
than  outside.  The  shear  at  the  boundary  has  increased  from  zero  for  the 
homogeneous  case  to  50%  of  the  normal  stress,  and,  as  seen  in  Fig.  18a,  it  is 
concentrated  outside  the  center  crystal. 

For  the  case  of  Model  B,  a  soft  Isotropic  inclusion  is  enclosed  in  a 
system  ten  times  stiffer.  The  stresses  on  each  face  are  reduced  by  a  factor 
of  0.27.  The  surrounding  structure  could  sustain  a  hole,  and  so  this  re¬ 
duction  is  easily  understood.  The  principal  compressional  stresses  inside 
the  center  crystal  are  less  than  0.20  times  the  compressional  stresses  out¬ 
side.  The  shear  is  raised  to  75%  of  the  average  normal  stress,  and  from 
Fig.  18b  it  can  be  seen  to  be  confined  to  the  region  outside  the  center  crystal. 
The  shear  inside  the  center  crystal  is  essentially  zero.  The  boundary  out¬ 
side  the  center  crystal  is  loaded  predominately  circumferentially. 

The  other  four  models  are  anisotropic.  Because  of  the  relationship 

v  E  =v  E  for  orthotropic  materials,  a  limit  is  imposed  on  the  degree  of 
xy  x  yx  y 

anisotropy.  The  ratio  of  3/1  permits  reasonable  values  of  Poisson's  ratio. 


Model  C  was  given  technical  constants  so  that  in  the  y-direction  the  constants 
are  continuous.  In  the  x-direction  the  modulus  is  three  tines  greater  than 
the  matrix.  It  might  be  expected  that  the  stresses  would  be  increased  on 
faces  one  and  four  (faces  perpendicular  to  the  x-axis)  and  that  the  other 
faces  would  be  affected  relatively  little.  The  finite  element  analysis  shows 
that  a  stress  increase  of  1.32  occurs  for  the  x-faces  an4  1.20  for  the  other 
faces.  The  geometric  effect  of  the  crystal  shape  on  the  trajectories  is 
evident  in  Fig.  18c,  and  it  is  evident  that  T  has  its  largest  values 
outside  of  and  all  the  way  around  the  center  crystal. 

Model  D  has  technical  constants  such  that  only  the  y-direction  is  different 
from  the  matrix  and  the  modulus  in  that  direction  is  1/3  that  of  the  matrix. 

It  might  be  anticipated  that  there  would  be  no  alteration  in  stress  at  the 
x-faces  and  a  reduction  on  the  other  faces.  The  analysis  shows  that  there 
is  a  stress  increase  on  the  x-faces  of  1.16  and  for  the  other  faces  a  re¬ 
duction  of  0.93.  In  Fig.  18d,  the  plot  of  the  principal  stresses  at  the 
center  of  each  finite  element  triangle  shows  an  unusual  behavior  near  the 
x-faces.  The  principal-stress  axes  are  rotated  nearly  90  degrees.  The 
values  for  the  principal  stresses  are  lower  inside  the  center  crystal,  but 
Tmax  *s  nearly  uniformly  large  inside  the  center  crystal.  The  direction  of 
principal-stress  axes  coincides  with  that  of  the  principal  strain  axes  out¬ 
side  the  center  crystal,  but  inside  they  do  not. 

There  is  very  little  difference  between  Model  A  (hard  isotropic)  and 
Model  E  (hard  anisotropic).  Figures  18a  and  18c  are  essentially  identical. 

The  stresses  extrapolated  to  the  boundary  differ  by  about  2%. 

It  appears  that  the  large  change  (xlO)  in  modulus  is  the  dominating 
factor.  The  stress-increase  noted  for  the  anisotropic  case,  where  the  inner 
elastic  modulus  was  three  times  greater  than  the  matrix,  was  1.32  (Model  C). 

For  Model  A,  where  the  modulus  change  was  ten  times,  the  stress  increase 
was  1.36.  Therefore,  the  three-fold  increase  in  the  elastic  constant  is 
seen  to  be  almost  as  effective  as  the  ten-fold  increase.  Hence,  the  aniso¬ 
tropic  effect  of  the  center  crystal  is  masked  to  a  large  extent  by  the 
larger  change  of  both  moduli. 

There  is  also  little  to  distinguish  Model  B  (soft  isotropic)  and  Model 
F  (soft  anisotropic),  by  examination  of  Figs.  18b  and  18f.  The  stresses  at 


the  boundary  are,  lioxvever,  reduced  on  all  faces  for  Model  F  and  on  the  faces 
perpendicular  to  the  x-axis,  the  forces  were  reduced  to  less  than  half  the 
value  they  were  subjected  to  for  Model  D. 

Six  Random  Orientation  Models 

In  Fig.  19  six  different  random  orientations  under  thermal  loading  are 
given.  The  center  crystal  is  maintained  constant  in  all  of  the  figures 
with  the  c-axis  in  the  direction  of  the  x-axis  and  the  a-axis  in  the  y- 
direction. 

The  position  of  the  maximum  stresses  varies  from  one  different  random 
arrangement  to  the  next.  There  is  a  tendency  for  maximum  stresses  to  develop 
near  the  corners  of  the  crystals.  The  stresses  extrapolated  to  the  boundaries 
from  the  inside  did  not  compare  well  with  the  extrapolation  from  the  outside. 
In  some  instances  the  extrapolations  from  one  side  of  the  boundary  were 
large  and  negative  and  from  the  other  side  were  large  and  positive.  A  finer- 
mesh  model  seems  indicated. 

CONCLUSIONS 

The  elastic  behavior  of  an  alumina  ceramic  body  has  been  modeled  by 
the  use  of  the  finite  element  method.  It  is  found  that: 

1)  The  anisotropic  characteristics  of  the  material  alter  the 
position  in  the  model  where  the  magnitude  of  the  shear  stress  is 
maximum  under  isostatic,  uniaxial,  and  thermal  loadings. 

2)  Under  uniaxial  loading,  the  effect  of  anisotropy  under  the 
conditions  examined  was  to  alter  the  stresses  acting  on  the  faces  per¬ 
pendicular  to  the  loading  forces  from  compression  to  tension  by  about 
6%  of  the  applied  loading. 

3)  The  inclusion  of  a  crystal  with  elastic  constants  ten  times 
greater  or  less  than  the  matrix  makes  such  a  difference  that  the  effect 
of  a  three  fold  anisotropy  is  nearly  masked. 

A)  Stress  gradients  generated  in  a  model  with  random  orientation 
of  the  crystals  is  very  high.  Further  conclusions  about  the  effect  of 
random  orientation  would  require  a  model  better  adapted  to  the  high  - 
stress  gradient. 
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Nineteen  crystal  model  of  a  two  dimensional  ceramic  body.  Each 
hexagon  can  be  given  independent  anisotropic  elastic  and  thermal 
constants.  The  orientation  of  the  axes  are  shown.  The  faces  of 
the  center  crystal  are  numbered. 


Finluo  element:  model  Cor  nineteen  cryotal  model  ofi  q  commie 
body.  The  model  lias  366  triangles  and  199  nodes, 
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The  elastic  behavior  of  a  ceramic  body  has  been  modeled  by  use  of  the 
finite  clement  method.  The  model  simulates  a  set  of  nineteen  hexagonal 
crystals  whose  elastic  and  thermal  properties  can  be  independently  varied. 
Several  different  configurations  of  the  elastic  properties  of  the  crystals 
were  evaluated  under  conditions  of  isostatic,  uniaxial  and  thermal  Icadings. 
The  effect  of  anisotropy  shifts  the  position  where  the  maximum  shear  acts 
in  the  model,  under  isostatic  and  thermal  loadings. 
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